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The quantum statistical treatment of the Rutherford model, considering matter as a system of 
point charges (electrons and nuclei) is analyzed. First, in the historical context, the solutions of 
ri . different fundamental problems, such as the divergence of the partition function, elaborated by 

Oj' Herzfeld, Planck, Brillouin and Rompe - most of the relevant papers published in the Annalen 

' ' der Physik, are discussed. Beyond this, the modern state of art is presented and new results are 

given which explain, why bound states according to a discrete part of the spectra occur only in 
rS , a valley in the temperature-density plane. Based on the actual state of the quantum statistics 

j^ ■ of Coulomb systems, virial expansions within the canonical ensemble and the grand ensemble 

'~^ ' and combinations are derived. The following transitions along isotherms are studied: (i) the 

Mh| formation of bound states occurring by increasing the density from low to moderate values, 

C/3 , (ii) the disappearance of bound state effects at higher densities due to medium effects. Within 

O ' the physical picture we calculate isotherms of pressure for Hydrogen in a broad density region 

' ly-j ' and show that in the region between 20 000 K and 100 000 K and particle densities below 10^^ 

^►s,, cm~'^ the cross-over from full to partial ionization may be well described by the contributions 

ri . of extended ring diagrams and ladder diagrams. 
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oo ■ 1 Introduction 

o. 

C*^ , Hundred years ago, Rutherford invented a new model for the interpretation of existing scattering 

experiments on the scattering of electrons on matter. In May 1911, Rutherford came forth with 

CO . a model for the structure of atoms which provided a first understanding why electron scattering 

^^ ' can so deeply penetrate into the interior of atoms, so far unexpected experimental results [T]. 

Rutherford explained his scattering results as the passage of a high speed electron through an 

atom having a positive central charge -|-iVe, and surrounded by a compensating charge of N 

electrons [l]. In Rutherford's model the atom is nearly empty, it is made up of a central charge 

and electrons and the glue keeping the system together are the Coulomb forces. Today we denote 

^ , the positive central charge as the atomic nucleus, though Rutherford did not use the term " nucleus" 

^_' in his paper. The central point charge is surround by a cloud of point electrons orbiting around 

the nucleus. The first quantum-mechanical treatment of the Rutherford model was given in 1913 

by Bohr. According to the Bohr theory the energies and radii of the orbits are (in Gaussian units) 

^e^ 2 ^^ mem+ 

2fi s^ fie TOe + m+ 

Here and in the following we denote the so-called main quantum number by s in order not to be 
mixed with the density n. According to the Bohr model there are infinitely many levels close to 



Corresponding author: c-niail: ebeling@physik.hu-berlin.de 
" e-mail: wolf-dietrich.l<raeft@uni-rostocl<.de 
" e-mail: gerd.roepke(§uni-rostocl<.de 



Copyright line will be provided by the publisher 



2Ebeling, Kraeft, and Ropke: On the quantum statistics of bound states within the Rutherford model of matter 

the series limit s — > oo. A first statistical theory of the Bohr model was developed already in 1913 
by the Austrian physicist Karl Herzfeld P] . Herzfeld detected that the Bohr model had a serious 
deficiency, the internal atomic partition function 



a(T)=^exp(-/3i?,):=^ 



s exp 



/ 



s^ksT 



I = 



lie 
2nV 



(2) 



is divergent (/ = —Ei. ionization energy). 

We see that the contributions to the internal atomic partition function for Hydrogen according 
to the definition by Eq. ^ increase as s^ (see Fig. [T]) and the sum can get any value in dependence 
on the maximal number s,„ and is divergent. This problem exists for all elements. It is inherent 
in any Coulomb system, e.g. the electron - hole plasma in excited semiconductors that will not 
be considered here. Here we deal mainly with Hydrogen which is the most abundant element in 
the Universe, its physical properties raised, in the past, many works. 
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Fig. 1 Left panel: The terms in the atomic partition function defined by Eq. ([2]) for two temperatures 
corresponding to /?/ = 5 (red and green curves above) and /?/ = 10 (blue and purple curves below). 
Up to terms of order Es^ — ksT the terms in the usual partition function decrease monotonically and 
then they start to increase as s^ leading to the Herzfeld divergence problem. For comparison we show 
the Brillouin-Planck-Larkin renormalized partition function with the terms of a convergent series. Right 
panel: The region in the density-temperature plane of Hydrogen where the formation of bound states is 
expected, see text. 



The development of the quantum statistics of plasmas is intimately connected with the re- 
quirements of a new scientific discipline - astrophysics. In order to bring some systematics into 
the existing observations, to understand why spectra appeared under definite conditions, a the- 
ory was required which provides the abundancics of atomic energy levels in stars. Such a theory 
was created on the basis of the Rutherford-Bohr model by the work of Eggcrt and Saha around 
1920. Their work is an application of thermodynamics to the Rutherford-Bohr model, which 
was inspired by the papers of the Nernst school. The ionization problem was studied within the 
chemical picture and using Sackur's and Tetrode's statistical expression for the entropy of gases. 
Eggert published a first approach in 1919 0], his paper came to the attention of the young bengalic 
physicist Saha, who succeeded to improve the approach considerably [31[S]. Saha considers, in the 
paper from 1921, Calcium plasmas as an example. Formulating his equation for Hydrogen we get 
for the fraction of ionized atoms a the relations 

^"" " " A3exp(/3/); A^ ^ ^\ -.. (3) 



nK- 



K 



y/27rfj,kBT 



The solution of this quadratic equations is given by 



a = ^(7), 7 = nK{T), 



A{x) 



tJ^ 



Ax-l] 



(4) 
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Considering an ideal mixture of electrons, ions and atoms we obtain for the pressure 

/3p = (1 + A{-i))n. (5) 

Wc call the function ^(7), which monotonically decreases with the density, the Saha function. 
At high densities where the mean distance of protons approaches the Bohr radius, this leads to 
an enormous pressure acting on the neutrals which will finally be destroyed; we note however that 
this limit is not described in the present approach. The region were bound states may exist is 
shown in Fig. [TJ 

Eggert and Saha used the thermodynamic approach of Sackur and Tetrode which needs only 
the binding energy, avoiding this way the problem of divergence of the partition function. For 
more accurate treatment of astrophysical tasks the most important problem was the divergence 
of the partition function of the Bohr atom. As already said above, the terms in the atomic 
partition function for Hydrogen according to the definition by Eq. ([2]) diverge as s^ (see Fig. [1]). 
The easiest way is to cut the sum at the lowest terms corresponding to \Es\ — ksT. This way 
of solving the divergence problem is due to the classical paper of Planck published in 1924 [6], 
later in 1938 applied to special plasma problems by Riewe and Rompe [7]. Developing Plancks 
approach Brillouin proposed in 1932 a more smooth way of removing the divergence which led to 
the formula 

00 

^bpl(t) = J2 A'^MPih^) - 1 - m/s^)]- (6) 

s=l 

This procedure was justified only in 1960 by Larkin [8] and subsequently by the present authors in 
collaboration with D. Kremp by using strict quantum-statistical methods [9l[ini[lIl[l2l[T3l[T4]. For 
later purpose it is interesting to note, that this so-called Brillouin-Planck-Larkin (BPL) partition 
function may be expressed (just after developing the exponential function) by a different series 
which contains the Riemann ^ function 

fc=2 

In the present paper we will give a systematic quantum-statistical statistical theory of matter 
within the Rutherford model, including in particular the equation of state. We concentrate on the 
so-called physical picture which stays within the Rutherford picture and avoids the introduction of 
chemical species like atoms, molecules etc. We note that in some earlier work the Saha theory has 
been proven to become asymptotically exact in zero-temperature and zero-density limit [15[ 1161 I17j . 
Here we restrict ourselves mainly to a region of practical interest, i.e., to temperatures between 
20 000 K and 100 000 K and particle densities between 10^^ and 10^^ per cm^^. In some of our 
earlier work the transition to a chemical picture was made following the principle of equivalence 
that bound states are to be treated on the same footing as free particles [TOl [131 HI] • This approach 
was quite successful in the description of partial ionization [131 [HI HOI [211 [22l [23] • On the other 
hand, one has to accept that in the region of high densities several serious difficulties appear 
which are mainly connected with the impossibility of a clear distinction between free and bound 
states in a dense system |24] . This is the reason why we stay here within the physical picture 
taking into account the important higher order terms [25| . In a similar spirit, but using different 
techniques, a consistent treatment of bound state contributions within the physical picture has 
been constructed already by Alastuey et al. [26l [27l [28] . We show here first a treatment of the 
Coulomb singularities in the framework of the virial expansions within the canonical ensemble, 
then we discuss the fugacity expansions and identify the higher order contributions which allow 
a Saha-type approximation. In the last part we investigate the role of the identity of electrons 
within the grand canonical ensemble. 
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In the first part of this work we will to show that the Brillouin-Planck-Larkin procedure provides 
the most natural way to avoid the divergencies for non-degenerate systems. In the last part we 
study degenerate systems and show that a consequent treatment of the identity of electrons in 
the Hartree-Foek approximation including bound states is a key for the treatment of hydrogenic 
bound states. The present approach is based on the method of Green's functions as developed in 

[Ml [HI mi Eg [30]. 

Any successful description of plasmas has to go beyond the Saha theory which is something like 
the lowest density order approximation for plasmas. This way our strategy should be to obtain 
first an approximation equivalent to the Saha description and then to go further. The derivation 
of Saha-type approximations from quantum statistics was first studied by Planck, Fowler and 
Brillouin [6l 1311 I32j . Since then it is a topic of permanent interest due to the importance of Saha- 
type equations for experimental work and for many technological applications. A first big progress 
in the strict quantum statistical treatment was achieved in the 60th by papers of Montroll, Ward, 
Vedenov, Larkin, Abrikosov and others [33l [34] and in subsequent work P [TO l [T2l [TS] [M l [T8 l 
[35l [36l [37l |38] . We note that a consistent description of bound states in macroscopic physical 
systems is a hot problem of modern quantum statistics and was treated also independently and 
using different techniques by other workers [26l [27l [28] . 

2 The model of Planck for Hydrogen plasmas 

Planck's work was inspired by the earlier papers of Eggert and Saha that was an application 
of thermodynamics and Sackur's and Tetrode's expression for the entropy. Accordingly Eggert 
and Saha needed only the binding energy, avoiding this way the problem of divergence of the 
partition function. In his Annalen paper from 1924, Max Planck uses a more consequent statistical 
approach based on the partition function. He splits first the atomic partition function into three 
contributions 

a{T) = a3{T) + a2{T)+aiiT). (8) 

where (JsiT) is the contribution of the lowest bound states \Es\ > ksT, (J2{T) is the contribution 
of the bound states below the series limit and cri(T) is the entire remainder. Using a rather 
complicated quasiclassical derivation Planck was able to show that due to compensation effects it 
holds 

a2{T)+a^{T)c^^. (9) 

In other words the sum of the second and the third contributions reduces to the partition function 
of an ideal gas, i.e. the all diverging contributions cancel each other. Planck's paper was considered 
in those years to be very important, so it was discussed in some detail in several textbooks and 
monographs, e.g. in the books written by Brillouin [32] and Fowler [31]. We follow here mostly 
Brillouin. According to Brillouin [32] , the cancellation effect found by Planck is exactly true only 
if one makes the particular choice 

^3(r) = aBPL(T). (10) 

As the derivations of Planck, Fowler and Brillouin are very difficult to understand we give the 
proof of Planck's important statement in a more easy version which however uses the same physical 
assumptions 0. We start with Planck's final expression 

F = kBTN[a[\n{nAla) + 1] + a[ln(nA^a) + 1] 
+ (1 - a)[ln(nA3(l - a)) - hi(a3(T)) + 1]] + F,^. (11) 
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Here Fcx is the excess part of the free energy. Planck and Brillouin claim that this contribution is 
zero, Fox = 0; however these researchers did not include screening, an effect detected nearly at the 
same time by Debye. In the following part of this section we repeat the derivation of the excess 
part but include also the contributions from screening effects using a semiclassical screened cluster 
expansion approach [9l |39] . In order to find the excess free energy of a plasma we have to calculate 
all contributions beyond the the ideal contributions of electrons, protons and atoms. Following 
[21 llOj and using the classical technique of screened virial expansions j39] we find including the 
Debye limiting law and the second virial coefhcents 



F,, = -knTV 



Bah = 27r i drr 



— — + 71gi?ee + ineUiBei + TL^ Bu 

LZtt 



Sab CXp(5a6 + fiVab) - I - gab - -^glb 



The Boltzmann factors without bound state contributions are given by [3] 

„2 /„^ T^ . _ „ „, /„ g^^ ^ Vabexp{-Kr) 



Sre{r) 



See{r) = Sii = exp(-/3e /r), Vab = eaeb/r, 

47r rvo 



^^^^ r ' V27r/ifci3T 



dpp' [eM-l3Eip, r))-l + l3E{p, r)] , 



Eip,r) = f-'-; 



Po 



2^ie^ 



K = 87rn/3e . 



(12) 
(13) 

(14) 
(15) 

(16) 



Here the second contribution in Sie is just subtracting the bound state contribution from the full 
Boltzmann factor in a classical approximation. Up to the Debye contribution, which was not yet 
known to Planck, since is was just in print when Planck wrote his paper, this is just what Planck 
and Brillouin calculated. After tedious calculations one can show analytically with a numerical 
check [9^, that the sum of the integrals beyond the Debye term cancel each other and gives zero. 



F,. 



-keTV 



K 

12^ 



+ 0(n'V2) 



(17) 



The excess free energy consists in the given quasiclassical approximation only of the Debye term 
and some higher order corrections. In particular, there is no contribution in the quadratic order. 
This way wc could justify the remarkable result of Planck and Brillouin, wc just found a correction, 
the Debye term. Including this essentially classical term we get 



/3p^{l + a)n - i(8^/3e2)3/2(an)3/2, 



(18) 



a = A{'-f) according Eq. (|3]). A typical picture of the pressure in the Saha approach shows that 
bound states exist only in certain valley in the density-temperature plane, see Fig. [1] We will see 
later that the bound state contributions as well as the Debye contributions are responsible for a 
reduction of the pressure in comparison to the ideal pressure (see Fig. [2]). The Saha equation and 
the corresponding pressure do not describe the region of high densities, the Debye term yields just 
a first correction. At higher densities another important influence is due to the Fermi pressure in 
the region of electron degeneracy. In order to check for the influence of electron degeneracy we 
compare with the ideal gas pressure including the Fermi pressure of the electrons: 



/3pf (/3,n) = n[l + 0.08839nA^-0.00083n2A6 + 0.000012n3A^ 



]; A, 



h 



\/2TrmJtBT 



(19) 



Looking at typical isotherms (see e.g. Figs. 2, 3) we see that the pressure related to the ideal 
pressure decreases with growing density and the pressure increases again. This demonstrates that 
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the region of bound state formation where the pressure is significantly smaller than the ideal 
pressure of an electron-proton gas is a limited region in the density-temperature plane (the valley 
of bound states). 

The main task of this work is the study of the two transitions: 
(i) the formation of bound states occurring at the cross-over from fully ionization to the region 
were bound states are formed (we have partial ionization) observed with increasing density, 
(ii) the destruction of bound states beyond the valley of bound states in the high density region. 

The region of our interest in the density-temperature plane where bound states appear is 
demonstrated in Fig. [TJ We determine here the first transition by the densities where p/2nkBT 
decreases below the value 3/4 and identify the second transition by the condition that p^ /2nkBT 
exceeds the value 1, that means the Fermi pressure p^ gives the overall dominant contribution. 
From the physical point of view we should expect that the real pressure should be nearer to 
the highest contribution, i.e at lower densities nearer to the pressure given by the Saha equation 
and at high densities nearer to the Fermi pressure. Later we will confirm the result that bound 
states occur only on an island in the T — n plane. In our primitive extended Saha model, the 
formation of atoms is described by the mass action law which is part of the Saha model and the 
destruction at high density is described in a very rough approximation only by the strong increase 
of the ideal Fermi pressure. In the following sections we will show how these effects are more 
correctly described by methods of quantum statistics. We present in Fig. [T] first a horizontal line 
corresponding to the density of n = 10^^ cm~^. For such a density the the mean distance between 
protons (or electrons) defined by 



47r 



-nr, 



1 (20) 



3 

is around the Bohr diameter. 

ro~2aB; na|~ 0.03. (21) 

For larger densities (smaller average distances) the formation of bound states may be excluded, 
the plasma behaves nearly as an ideal Fermi gas since there is not enough space for forming atomic 
orbitals or molecular orbitals. The disappearance of bound states is referred to as Mott condition, 
corresponding to the fact that atoms are destroyed if the mean distance of the electrons is much 
smaller than the Bohr radius. Such situation was discussed in [501 HOI EI]- For more detailed 
information see [Ml [3^. We show, in Fig. [H a line a ~ 0.5. The Saha equation provides for the 
condition a = 0.5 the estimate 

70 -: nK{T) = nA^ eMl/ksT) = 2. (22) 

Wc arc here mainly interested in the region of atomic bound states. Near to the upper border 
of the corridor region shown in Fig. [1] only the atomic ground state exists, but near to the lower 
border many bound states are to be expected and we have to avoid the divergence of the partition 
function by appropriate renormalization. At lower temperatures T < 30 000 K we may expect 
also the formation of molecules in a certain density region. 

3 Derivation of the Brillouin-Planck-Larkin partition function by 
low-density virial expansions 

Following our earlier papers [121 ES] a-nd books [13l [14] we use the method of cluster expansions 
in order to derive the exact virial functions. We start here with a treatment of the Coulomb 
contributions in the framework of the virial expansions within the canonical ensemble. For a real 
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quantum gas with only short range forces the free energy may be at low density described by a 
virial expansion [T^l 155] 



F = F; 



id 



keTV 



y^ naUbBab + ^ UaTlbncBabc 



ah 



Bab = const Tr[cxp(-/3i/,b) - cxp(-/3H^b)], 
where the Hamiltonian of pairs of species a and b is defined as 



(23) 
(24) 



H„ 



'^fJ-ab 



A + Vab. 



(25) 



Using the resolvent representation for the exponential operator we get 



Bab 

F{z) 



{2Sa + \){2sb + l)2TTi 
1 1 



exp(-^z)F(z)(iz, Xab = 



\/2^LabkBT' 



Tr 



Hab - Z i/ ' - Z 



(26) 



After some transformation we arrive at the following representation by Jost functions Di{z) 



F(.).(^±i±ii(|fi±iif:(2.+i) 



(1 + 5ab) 



e=o 



l±Sab 



hiy 



(2s, + 1) 



(-l)-lni^.(.). 



(27) 



The Jost functions Di{z) are analytical functions with poles at the bound states and a branch 
cut at the positive real axis defined by the scattering phase shifts. These functions as well as 
other scattering quantities are exactly known for Coulomb systems. However in the procedure of 
calculating the integrals there is a difficulty due to the long range character of Coulomb forces and 
the necessity to introduce some screening procedure. In analogy to the classical case discussed in 
the previous section we may write including screening effects 



F-F,'! 



keTV 



^3 _ 

— + ^ UaUbBab + ^ naUbricBab 



abc 



(28) 



Here the tilde denotes that screened potentials were included into the virial coefficients. In order to 
reduce the difficulties connected with a complete derivation we proceed as follows: The screening 
procedure is necessary only for the lowest orders of the interaction parameter e^. Higher orders 
in e^ beginning with e^ have already a decay rate of e^/r^ and do not need any screening. Taking 
into account that screening refers to the low orders we split the second virial coefficients into two 
parts 



Bab — G Bab + G Bab 



G Bab — G Bab 



(29) 



G'Bab = 0(1) + O(e') + Oie*) + 0(6^), G"Bab = 0(e«) + 0{e'°) + 



(30) 



Let us first calculate the higher orders, remembering that according to Brillouin and Larkin, 
here the contribution of the bound states are to be expected. Introducing here the known Jost 
functions, the final result without symmetry effects reads 



G"Bab = 27r3/2A3jl + peaebu] J2 



cm /i __ o2 — Tn\ 

-^ AQim - 2) ± 5ab- -C{m - 1)1 . (31) 

2'"r(m/2 + l)^^^ ' "'' 2sa + l Jl \ I 
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The contributions in lower order have to be calculated term by term. The result of tedious 
calculations is (C = 0.577 - Eulcr's constant) 

+ I (/3eaeb)3 [(1 + Peae^n) ln(3KA„b - /3eaebK(l - ln(4/3)] 



±5ab[l + fieaCbK] 



V^ jab V^(ln2) 2 ,^.-6 



f3 

72' 



(32) 



Note that most difficulties are connected here with the term of order ^^^ [T31 [HI [31] . 

For charge and mass symmetrical systems (all masses are equal and the reduced masses are just 
half of it) the odd contributions cancel and we get in the case of Hydrogen for the sum, taking 
into account (31 = ^^/4: 

^ G"B,, = 2nV^Y. ^^'^^^klf = 47rV^AiaBPL(r). (33) 

ab k=2 

This way we have given a strict derivation of the Brillouin-Planck-Larkin partition function given 
by Eq. (0 . The final results for the equation of state of Hydrogen may be expressed in form of a 
density expansion : 

- 2n2[AVBPL(T)][l + ^/3e2K)]-2n2iC*(T)]+L(T)n5/2 + 0(7i3logn). (34) 

The function K* (T) which appears in the formula Eq. ([34)) depends only of the temperature. This 
transcendent function is exactly known |13[ 1141 136] and may be calculated as an infinite series in 
^ using the expressions for Bah given above. There are also tables for Hydrogen available [13]. 
However since the full expression is a quite long and not easy to handle, we will give here only a 
rather good (in numerical respect) and simple approximation 

/r(r) = -iA^/3/-^A3. (35) 

This expression is exact with respect to the largest asymptotic term 0{I /ksT) and was checked 
with the available tables [13]. As a technical remark we note, that the terms not contained 
in the approximation eg. ([55]) are mainly given by the contributions of the odd powers in ^ to 
the pressure. We do not give here a detailed discussion since these terms give only a small 
contribution to the pressure in the region lO^K < T < 10^ K. The exact expression may be found 
in original articles [TUl [H] and in several books [T31 [HI [3S] . The complete term of order rv'/'^ 
has also been computed in several works [37l [42l [43] . We will not discuss this matter here since 
we use here some approximations for L{T) and because this function does not contain bound 
state contributions. There were different attempts to determine such contributions n^/^ in the 
framework of Green's function techniques. Here we mention jST] [38]. These results should be 
compared to the corresponding expressions for the n^/^ terms given by [331 [H] ■ Though there are 
results which seem to be in agreement with each other, there remains still some open questions 
[4T| , see [45] . Here we are interested mainly in the bound state effects and restrict ourselves to the 
asymptotic results for the case ^^^ » 1 what is equivalent to T ^ I /kB- For the case of Hydrogen 
and lower temperatures T <^ I/ks we get in asymptotic approximation the final results for the 
pressure at small densities: 

8v 2 2 

- n^A'^{l3I)+n^A^aBPLiT)[l + ^l3e^K] + --- + Oin^/^) + 0{nHogn). (36) 

Copyright line will be provided by the publisher 



andp header will be provided by the publisher 



This way we have shown that in the quantum-statistical expressions for the thermodynamic func- 
tions appears the Brillouin-Planck-Larkin partition function in a natural way as a convergent 
infinite sum starting with terms of order e®. As mentioned in the previous section, a first approx- 
imation to this partition function was obtained by Planck [6] and a first explicit derivation was 
given by Brillouin j32j . Later an asymptotic quantum-statistical derivation was given by Larkin 
[8]. We follow here Larkin's asymptotic approach which was extended later in [10] but include 
more terms than in the earlier work. Exact expressions for the orders 0{n) and ©(n^/^) in the 
pressure including all orders in e^ were derived by using cluster expansion method and Green's 
function methods in [10l[l2l[l3l[l4l[35l[36] and by using different techniques in [26l[27l[28]. We will 
not give a complete survey about the existing exact results about the lower orders in the density 
but will instead direct the following study on the most important higher order terms. However 
let us still make a remark about the foundation of the Brillouin-Planck-Larkin partition function. 
There seems to be a never-ending discussion about the Brillouin-Planck-Larkin partition function. 
Just to give one example, Starostin and Roerich [46] derived recently a quite different expression. 

Our point of view which is based in particular on the results obtained here is the following: 
(i) The quantum-statistical result for the contribution 0{n'^) to the pressure is exact. The result 
of several groups working with independent methods give the same results (see e.g. [Ml [26]). 
(ii) The BPL-partition function is asymptotically in full agreement with the second virial co- 
efficient. There is a small freedom in the special choice of a partition function, as far as the 
asymptotics is not violated. However so far no good reason is to be seen, why we should leave the 
BPL- partition functions and make another choice ^^. 

(iii) According to Onsager's bookkeeping theorem, in a correct theory any small change of the 
partition function does not matter, the final results for any physical observables are stationary 
with respect to small changes [24] . 

Some of the critics of the BPL partition function is just based on intuitive arguments and is not 
based on a serious revision of the quantum-statistical foundations given in the work of Vedenov 
and Larkin and the subsequent papers. 

Admittedly the quantum statistics behind this formula is very complicated but nevertheless it is 
strict. With respect to the intuitive arguments we only want to say, the formal reasons for having 
no terms of order 0(e°) and 0{e'^) in the Planck-Brillouin-Larkin formula is that these terms were 
consumed already in the formulation of the ideal terms which are 0(e°) and the screening terms 
0{e^), and C(e*) which arise formally by summing up an infinite series in ©(e"). This procedure 
consumes some terms of order 0{e^), and all contributions of order ©(e*) to the partition function 
which cannot appear anymore in exact expression for the partition function, if double counting 
of terms is strictly avoided. Another, formally equivalent argumentation is based on the fact of 
compensation of the terms below the series limit and above the series limit is most essential in this 
context [36] . We want to underline again that there is no unique choice of the partition function 
for Hydrogen, there are other possibilities [24]. Essential is however, that all possible choices 
should be compatible with the quantum statistical formulae given above. That means among the 
many admittable versions of partition functions, the essential criterion for correctness is Onsager's 
bookkeeping rule. In other words the question is: Is there a correct counting of all contributions 
or not. Is there a check for double-counting or not [24]. As a matter of fact, all partition functions 
which work without a BPL- or equivalent subtraction of first order terms may have the problem 
of double-counting of certain diagrams of order 0{e"). 

Using the asymptotic expression for K*{T), which is asymptotically exact for ksT <C /, the 
pressure may be written in a different form 

/3p(/3,n) = „[l + _l=„A^ + ...]+n-^[l-^M) + l(^Af + ...] 

- n2[AVBPL(r)][l-/3e2Ac]+0(n3logn). (37) 
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Fig. 2 Pressure related to the total classical pressure 2nkBT at constant temperatures T = 0.1/ and 
T = 0.15/ in dependence on the density (log-scale). The result of density expansions up to 2nd order 
virial terms (red and green) are compared with corresponding fugacity expansions. 

We combined one of the two contributions from K*{T) with the ideal electron pressure and the 
other one with the Debye law. Further we introduced the mean thermal wave length A = A^e 
which is an effective quantum length, in some correspondence to the Debye-Hiickel parameter in 
the classical theory (/x - reduced mass). We have to note here again that the results given above 
are not fully complete in the order n'^/^ and in missing higher order corrections. However this 
problem seems to be not relevant for the bound state problem. A formally equivalent form of Eq. 
([57]) which shows better the general structure is the following 



/3p(/3,n) 



n[l 



8\/2 



nA^ H ] + n - ?17[(1 - Pe^n -\ ) + 



^[l-|7 + -l|l-2^(»A, + l(»A,» + 



■]+0{n^'^) 



(38) 



7 = AVBPL(r) = 87r3/2A3^f7BPL(r). 



(39) 



Summarizing wc may state that there exists an asymptotically exact expression for the pressure 
of Hydrogen in the limit of small densities, which is given in the formulae Eqs. ()37|) . p8|) which 
we have drawn in Fig. [2l This approach is compatible with the extended Saha theory what 
encourages us to proceed in this direction. Let us first investigate Eq. p8|) as it is, without 
completing the series. Looking at the curves in Fig. [2] we see that the screening effects and the 
corresponding BPL effects tend to lower the pressure. We see that this lowering is too large at 
increasing densities. 

The structure of Eq. (|38p suggests that the expressions in the parenthesis may be only the first 
terms of some infinite series with respect to nA^, k and 7. Just to give a first example: The first 
series is evidently nothing else than the ideal electron pressure, which may easily be extended to an 
infinite scries using the known expression for the ideal Fermi pressure. There is no question that 
all these terms appear in the full cluster series. In order to find the higher orders in the other series 
we have to use a more complicated procedure which is based on fugacity expansions |10l I13( I47j . 
According to our general philosophy, the strong decay of the pressure observed in second order 
density expansions (Fig. [2]) is due to the fact that important higher order terms corresponding to 
some infinite series were omitted in the derivation based on the canonical ensemble. To circumvent 
this difficulty is not an easy task. We mention that Alastuey et al. [27l [28] succeeded already 
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in deriving a Saha type expression in summing up an infinite scries of terms. Here we follow 
a similar idea but go an alternative way through fugacity expansions. As shown much earlier, 
fugacity-like expansions are equivalent to mass-action laws. This was used in plasma theory by 
Bartsch, Ebeling and others [13l [TOl |47] and worked out in large detail by Rogers et al. [48] . We 
mention that Rogers et al. succeeded to show that extended fugacity expansions allow an excellent 
description of measurements for the oscillation modes of the sun [49] . 

4 Combined density-fugacity expansions including nonlinear ring 
and ladder contributions 

At higher densities a systematic quantum statistical approach to the pressure of plasmas can be 
given based on the Green's function representations of the pressure. The central relation for the 
pressure reads [TS] [TU [18] 

p{l3,^le,f^^)=P^d-^ J — J dldlV {li)G2{lA,l + + , 1+ : h = 4) , (40) 

where 1 = {pi,(7i} denotes momentum and spin variables. The Green's function representation 
works in the grand canonical ensemble which provides us a series in the fugacities 

Za = —^3 CXp(/3^a) = 7la CXp{f3lll^) 

■''■a 

instead of a series in the densities. Note that fi"^ is the excess part of the chemical potential (the 
part beyond the ideal Boltzmann term) and that therefore Za — > n-a for small densities. At small 
fugacities (densities) the pressure is given as an expansion in the fugacities 



/3P = /3 ^ ^a + ^ ZaZbbabiHg) + ^ ZaZbZcbabc{Kg) + ^ ZaZbZcZdbabcd{Hg) H , (41) 

o. ab abc abed 

where 

K2=4^/3^z,e2 

a 

is the grand-canonical screening length. The structure of this series is similar to the density 
series, we have e.g. similar screening effects, however there are some differences, e.g. the fugacity 
series contains more diagrams. The contribution of bound states of an electron and an ion (atom) 
is contained in the contribution ZeZib^i, the contribution of molecules is contained in the term 

General expressions for screened fugacity series were written down first by MontroU and Ward, 
Vedenov and Larkin and explicit calculations were done by DeWitt and Larkin [S| 1331 [10] • Semi- 
classical fugacity series were given by several workers [13l [47] . 

We are interested here mostly in the bound state effect and may simplify our calculations 
by using approximations allowed for the region T < 157 000 K where the bound states appear. 
Further in order to simplify the quite complicate evaluation of these diagrams we restrict the 
theory to the region of non-degenerate plasmas. As to be seen from the above formulae, the 
difference of the masses of electrons and ions (protons) do not play a role in the region \£^ab\ ^ 1 
and the plasma behaves asymptotically like a system with equal masses and equal relative De 
Broglie wave lengths A = Aie- In the lowest approximation corresponding to Eq. (|37p we find for 
low temperature Hydrogen plasmas the following results valid up to order 0(z^/^) [T3 l ll4 j . 



^3 

I3p^ze + Zi + -77-/(%A) + 87rZeZ,;A^ exp[l + /3e^Kg]aBPL{T) + ■■■ , (42) 
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where the grand-canonical screening quantity is now defined by 

Kg = 47r^(ze + Zj)e^, 
and the quantum-statistical ring function is a transcendent function |10[ 113] with the series 



/(^) = 



1 3v^ 1 2 

1 X H X 

16 10 



(43) 



The relations between densities and fugacities are given by 



d(3p _ ^l d 



z,^r- = z, + -^i^iKgXfiKgX)) + Anz,z,[l + Pe^Kg -f • • • ]AVbpl(T) • • • , (44) 



dBv K^ d 

z,^ = ze + T|TTT(^gA/(KgA)) + 47rzezai + Pe" Kg + • • • ]AVbpl(T) • • • . (45) 



We will show that going by iterations from the fugacity variable to the density this result occurs 
to be equivalent to eq. (j37p . In order to proceed with this complicated system of equations we go 
to the special case of non-degenerate hydrogen at lower temperatures T < I /kg and satisfying 
nA^ < 1. As shown above, in this region the differences of the electron ion masses does not matter 
and we my assume z,, — Zi = z. In this approximation, Eqs. (j42p . (|44p . (j45p read 

I3p = 2z+ -f /(kA) + 8^z2aVbpl(T) + • • • , (46) 

" = ^ + TTT tt('*s^/(^9^)) + 8'^^'A-Vbpl(T) • • • . (47) 

iOTT OA 

In order to represent the pressure by densities, we neglect, in a first step, the bound state contri- 
butions. Expressing first the fugacities by densities we find 

z = nexp[-^G(A.A)], (48) 

where HUl H 

Gix) = l-^x + lx'----. (49) 

Then excluding the fugacities from the series step by step we arrive at the representation 

Pp = 2n-^^iK\). (50) 

Here 0(x) is another transcendent function, related to /(x) and /'(x), the so-called pressure ring 
function which is an analogue of the Debyc-Hiickel ring function [131 [33l [50] . The first terms of 
the series expansion read 

0(x) = [l-^x+lx^-...]. (51) 

A full representation of this quite complicated transcendent function may be obtained by methods 
similar as described in ]13j . We found the infinite series 

m = l-—^ ^y +3^ 2^(2fc-l)!! ■ ^^') 

fc=l k=l ^ ' 
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Fig. 3 The nonlinear ring function (^{x), which determines the contribution of screening to the pressure 
(blue line). Further we show the nonlinear ladder function A{x) (red line) which determines the contribu- 
tion of the bound states to the pressure as well as the linear approximation, which gives only the lowest 
order to the second virial coefficient (thin line). 



The asymptotic behavior of this function is given by 



H^) 



(53) 



This suggests the convenient Pade representations, which will be used in the calculations 
24 „_ 1 



0(x) 



G(x) 



l + (0F/4)x' 



24 + 90Fa; + 8x^ ' 
A representation by error functions reads 



(-) = ^ + ^ - ^ exp[xV4] - #[1 - <J>(x/2)] + ,F,. 



3x 



3.T 



(54) 



(55) 



The saturating behavior of the nonlinear function $(x) in comparison to its linear approximation 
is shown in Fig. [3l 

We summarize the results we obtained so far: By means of grand-canonical methods we suc- 
ceeded to obtain some infinite series in the screening parameter k instead of the first linear and 
quadratic terms only, as we found in the density series, (see Eq. p8[) ). This summing up a series 
in K which leads to a saturating function (see Fig. [3]) and improves very much the behavior at 
larger densities as we will show. 

In our representation Eq. (j38p . there appears also another series expansion in the density 
parameter 7 . In order to sum up this series in 7, we proceed in a similar way. In order to 
understand the structure of the series we study at first the bound state contributions, neglecting 
the contributions coming from the ring term and from the terms due to degeneracy. Further we 
replaced the factor 1 — /Se^n by a more smooth exponential function exp(— /3e^K). This yields the 
simple quadratic equation 



ttz'AVbplIT). 



(56) 



By the way, this formula shows already that the fugacities of the electrons and the protons should 
be rather small in the bound state region where the partition function is large. The quadratic 
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equation for the fugacities can be easily solved. This way we find, after elimination of the fugacity 
z (alternatively using the A - function): 

1 



zO = n^(7); A(x) = — [^/TTl^ - 1]; ^i ^ %T:n\^ g^^^{T) . (57) 

We note the following series expansion and large x asymptotic of the nonlinear function A(x): 

A(x) = 1 - a; + 2a;2 ; A(x) = ^ + • • • . (58) 

The saturating behavior of the nonlinear function A{x) in comparison to its linear approximation 
is shown in FiglSl This shows that in the region of large partition functions the fugacities disappear 
as 

0° = ^ (59) 

Y'87r7iA3crBPL(r) 

This is important for the understanding, why the fugacity series has, in the region of bound states, 
a better convergence than the density series. The better convergence is due to the fact that the 
fugacities disappear in the region where the bound state contribution is large. 
Introducing the zeroth step of iteration z" into the pressure we find 

/3/ = z" + (zO + 8^z2A3aBPL(T)) ^ z" + n = n{\ + ^1(7)). (60) 

This surprising representation tells us that the Saha function A(x) appears in the pressure in a 
quite natural way through the fugacity expansion. We have introduced elements of the chemical 
picture without using the notation of chemical species explicitely. The new " semi-chemical rep- 
resentation" includes the nonlinear function ^(7) which saturates at large densities (see Fig. [3]). 
The new representation contains all terms in the density up to r?' as well as several higher order 
terms in the density contained in the nonlinear ^(7) which is related to the Saha function defined 
by Eq. ([4]). Wc need this kind of functions to reproduce an ideal Saha-type behavior. The fugacity 
series contains more terms than the density series, so we may expect that some of the difficulties 
connected with density expansions, as the strong decrease of the pressure with increasing density 
shown in Fig. [5] may be avoided. Indeed, a representation of the curve corresponding to Eq. (|42p 
shows a more reasonable behavior with increasing density which is much closer to the Saha-type 
behavior. As the curves shown in Fig. [2] demonstrate, the pressure according to the fugacity ex- 
pansion goes to saturation. This is due to the fact that the fugacity is not fixed, it decreases with 
increasing value of the partition function and this way limits the growth. Evidently the fugacity 
expansions provide a more correct description of the bound state contributions as demonstrated 
in Fig. H 

We find that density as well as fugacity expansions have both advantages as well as disadvan- 
tages: 

1) The density expansion describes well the screening effects but it fails to cope with the diverging 
contributions from the screening terms and from the BPL-partition function. 

2) The fugacity expansion corresponds to an infinite density series including the partition func- 
tion <j. If a is large, then the fugacity goes to zero what guarantees even at large densities always 
finite contributions to the pressure. This is true for the screening contributions and for the bound 
state contributions; any strong increase of contributions suppresses the fugacities. In general, the 
fugacity expansions are more smooth. 

In conclusion wc may expect that the best representation is obtained by extended density 
expansions which contain additional contributions corresponding to the important damping terms 
in the fugacity expansions. We may expect that this procedure provides Saha-like terms. This 
kind of mixed expansions combines the positive features of both expansions avoiding the negative 
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features. In order to demonstrate this we started here from the density expansions and used the 
fugacity expansions mostly only for finding the right continuation of the infinite series with respect 
to the 7- parameter and the kX parameter. 

Following this line and including all terms up to known orders in the densities and fugacities, 
we obtain, in a first step, the following expression for the fugacities 

z = na'; a ^ A (STmA^ exp[-/3e2KG(KA)]crBPL(T)) . (61) 

Note that here the second term plays the role of an effective density of the bound particles. For 
the pressure we find, including again the full electronic Fermi pressure, the following relatively 
simple formula which contains the three nonlinear functions which we discussed in the foregoing 
part 

^p = pf + nA{V) - ^0(«'A), (62) 

with the definitions of renormalized values of the Saha parameter and the screening parameter 

r = 87rnA^ exp[-/3e^KG(K'A)]crBPL(r); k' = kA^/'^{-i). (63) 

This provides us with a closed and relatively simple formula for the pressure, including the Fermi 
pressure. The underlying assumption is that all electrons contribute to the Fermi pressure. This 
result for the pressure according to the combined density-fugacity representation Eq. (j62p which 
includes all contributions up to second order in the fugacity is, for three temperatures, shown in 

Fig. m 

The good convergence of the new expression including the nonlinear functions A{x) and 4i{x) 
based on the fugacity series is explained by the fact that the fugacities of the electrons and 
the protons are rather small in the bound state region. We notice again that this result comes 
from an extension of the density expansion including all quadratic terms from the fugacity virial 
expansions and the full Fermi pressure of the electrons. We see that the overall behavior of the 
new representation is much better than that of pure density or pure fugacity representations. This 
way we may conclude that the most appropriate description of Coulomb systems is by density 
expansions and including some elements from the fugacity expansions. The physics behind is 
that the fugacity expansions describe well the saturable forces between bound states in Coulomb 
systems, but on the other hand the additive long range Coulomb forces and the screening effects 
are better represented by density series. 

We note that the new theory based on Eqs. (p^ . (|55| is consistent with the Saha theory and 
in particular also with the Saha-Debye-Hiickel theory [13]. Comparing the curves based on the 
extended theory. Figs. |4] and [5] with Fig. [21 we see that the overall behavior of the new extended 
density representations is much better than that of pure density or pure fugacity representations. 
Further we may state that the new extended analytical theory is, up to the density n ~ 10^^ 
cni^^, in quite good agreement with earlier calculations based an the chemical picture and quite 
complex numerical codes [T9l [2T1 [22| [24l [51] 

Note however that the agreement is based on a complete agreement in the lowest orders in the 
density, in higher orders several differences appear. Beside the good agreement with the Saha 
theory and in particular with the extended Saha-Debye-Hiickel theory developed in J13j we list 
the following deviations: 

1. The exponent exp(/3/) in the Saha equation is replaced by the BPL-partition function in the 
quantum-statistical formula Eq. (|62p . 

2. The linear first Debye-Hiickel screening function linear in n, is replaced here by a nonlinear ring 
function (j){x) describing the quantum statistical ring sums in the grand canonical expansion in a 
good approximation. 

3. The linear second virial coefficient is replaced by a nonlinear ladder function stemming from 
the representation in the grand ensemble. 
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Fig. 4 Equation of state (EOS) calculated from an extended density representation including the non- 
linear ring and ladder functions (\)(x),A(x) and including the Fermi pressure of the electrons. We show the 
pressure related to the classical ideal pressure inksT at constant temperatures T = 20000, 30000, 50000 K 
in dependence on the density (in log scale). For the temperature T = 20000 K (lowest, green line) we give 
a comparison with an earlier calculated curve based on the chemical picture and and Fade approximations 
Il9l[2Tl[22]. 

We note that the grand canonical ensemble plays an essential role in our derivation, and further 
we note that the formula Eq. (j62p corresponding to a second iteration is in some sense incomplete, 
e.g we would expect that in higher iterations more K-tcrms arc replaced and that more terms 
corresponding to a mass-action law will appear. 

We underline that several of the terms beyond r? and z^ are based on extrapolations which 
still need further confirmations. However in the region of low temperature T < I /ks and non- 
degenerate plasmas our rather simple formulae give a rather good behavior and describe well the 
transition from low density to the valley of bound states. We note however that several physical 
effects as e.g. plasma phase transitions [131 HI] are not yet described by the present approach. 
Evidently this effect appears only in higher order terms or after transition to some chemical picture 
p!3l [19] . However this is not our aim here. 

In order to describe the transition to full ionization in the degenerate region as well, additional 
effects have to be taken into account, in particular the symmetry between electrons and protons 
is lost. The disappearance of the bound states with increasing density was discussed in [^ HT) : 
for more information see [T31 [35] . 



5 Discussion and Conclusion 

It is shown in this work that transitions from full to partial ionization in the temperature range 
10^ — 10^ K which are due to the formation of H-atoms, are well described in the physical picture. 
Ring and ladder diagrams in higher order approximations as well as certain combinations of these 
basic diagrams are taken into account. In particular the divergence of the partition function of 
the Bohr atom and the state of art concerning this problem is discussed. 

Formally the divergence of all partition functions within the frame of the Rutherford model 
is due to the infinite range of the Coulomb potential. Practically, a divergent partition function 
means that ionization is impossible, atoms are absolutely stable. As stated already by pioneers 
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like Herzfeld, Fermi, Planck and others, this makes no sense and there should be physical reasons 

which make the atomic partition function finite. Among the physical effects leading to a cut of 

the partition function at some level, which were mentioned already by the pioneers, are: 

(i) The high levels \Es\ < kBT are not stable against collisions in a thermal system and are 

strongly perturbed by screening effects, this is represented by the BPL reduction. 

(ii) At higher densities and moderate temperatures, the contribution of the partition function 

provided by the ladder contributions strongly increases, also the nonideality strongly increases. 

The total contribution of screening and bound states remains however finite since the fugacities 

of the charges become small. The nonideality contributions saturate and lead to a reduction of 

the pressure to about one half of the ideal pressure (see Fig. 4). 

(iii) At very high densities beyond 10^'^ cm^'^ the electron degeneracy delocalizes the electrons and 

makes atoms obsolete. The nonideality decreases again, the plasma leaves the region of nonideality 

(screening and bound states). 

The principal scheme how the formation of bound states depends on the density at a fixed 
temperature is shown in Figs. [5] and [6l In Fig. [5] a comparison of the theory given here for 
T = 20000 K, T = 50000 K with the numerical results obtained within an advanced chemical 
picture by minimization of the free energy is presented. With the chemical picture we mean here 
a description which starts from an expression for the free energy depending on the density of free 
charges, atoms and molecules. One can see that the overall agreement is quite reasonable, the 
deviations increase only at large densities beyond n ~ 10^^ cm""^. The largest deviations appear 
near the minimum of the relative pressure for Hydrogen at 20000 K and may be interpreted 
by the formation of molecules which were not taken into account in the present work. These 
interpretations are supported by the the last Fig. [51 Here we show the degrees of ionization 
and the degree of molecular bound states obtained within the chemical picture for T = 20000 K. 
The composition and state of the system was obtained by minimization of the free energy. The 
numerical results used in Figs. [5] and [6] for comparison were obtained with a code developed in 
some earlier work (24l [30l |40] . 

Looking at the pressure curve in Fig. [5] one observes as in earlier sections first a decrease 
corresponding to the formation of bound states (atoms and molecules) . This decrease comes to 
saturation and we observe a valley of developed bound states. With increasing density the spatial 
occupation of the atoms and molecules as well as effects of electronic degeneracy come into play and 
the bound states are destroyed. The largest deviations between the curves for 20000 K are due to 
the formation of H2 molecules (deep minimum in the blue curve) which are not taken into account 
in the present analytical theory and also effects of strong degeneracy. These rather complicated 
effects are not discussed here. Parallel to the destruction of the bound states we observe an increase 
of the relative pressure in the region of rather high densities n > 10^^ cm~'^. The monotonic 
increase of the relative pressure may be interrupted at temperatures below 10000 K by a small 
wiggle. This is not an artifact. A rather fiat region of the pressure in dependence on the densities 
is contained. At lower temperatures the wiggle may be even stronger and should be discussed 
in connection with the possible existence of plasma phase transitions [13l [TH [20l [21] [22l [29l [51] . 
Recent work in this area was reported in |541 I55j and found a phase transition. This is not our 
topic here, the purpose of the present work is to consider the region of lower densities and higher 
temperatures. Our aim is, to obtain similar results as shown in Fig. [5] without referring to the 
rather complex chemical description based on atoms and molecules and their interactions. Instead 
of the many assumptions contained in any chemical description, the calculations given here are 
based only on Rutherfords model of point charges and Coulomb interactions as well as on the 
exact principles of quantum statistics. No a priori knowledge about the properties of atoms and 
molecules is required. The price to pay is that we have to avoid the regions where these species are 
dominant. With respect to accuracy the rather simple formulae given here may not compete with 
the recent successful approaches to dense plasmas based on numerical methods [511 [S31 [M] ■ We 
are convinced however, that in particular to experimentalists, analytical approaches may always 
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Fig. 5 Comparison of the analytical theory for T = 20000 K (magenta) and T = 50000 K (green) to 
the results of a numerical code developed earlier within the chemical picture [211 1221 1241 1511 I53j (blue: 
T = 20000 K, red: T = 50000 K). We show the relation of the pressure to the classical ideal pressure of 
the ions nksT. The largest deviations between the curves for 20000 K are due to the formation of H2 
molecules (deep minimum in the blue curve) which are not taken into account in the present analytical 
theory. 
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Fig. 6 Typical isotherms of the degree of ionization at dissociation of hydrogen for T = 20000A' in a 
wide density range based on the code developed earlier within the chemical picture [211 1221 1241 [51] . We 
show the degree of ionization (upper curve) and the degree of protons bound in H2 molecules (transient 
maximum) as a function of density in log-scale. 



be quite useful in complementing numerical results. In conclusion, we first analyzed the basic 
results for the quantum statistics of the Rutherford model of matter, partially in the historical 
context starting with the work of Bohr, Herzfeld, Saha, Planck, and Brillouin. Then wc discussed 
the essential results of quantum statistics since the 60th of the 20th century. The main new 
results obtained in this work are: For Hydrogen plasmas at lower density the EOS within the 
physical picture was obtained using a combined density fugacity representation which contains 
the essential contributions describing Saha effects and nonideality. The good convergence of the 
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fugacity series is explained by the fact that the fugacities of the electrons and the protons are 
rather small in the bound state region. The obtained new expression for the EOS is valid for 
non-degenerate plasmas and temperatures between 20 000 K and 100 000 K. Partial ionization, 
i.e., the appearance of atoms, is restricted to a valley in the density-temperature plane. There is 
a transition from below starting in the region of nondegenerate electrons and a transition from 
above beginning in the region of degenerate electrons. Both transitions are connected to different 
physical effects. 

I has been shown that in the non-degenerate region, the most essential effects is the BPL- 
reduction, resulting in the combined action of collisions and screening by absorbing the orders 
in e^,e'^ and the compensation effects of the states below and above the series limit. The most 
essential new contributions given here are the nonlinear ring and ladder functions introduced 
into the EOS. These terms stemming essentially from the grand canonical representations are 
responsible for a good representation of the EOS in a wide range of densities and temperatures 
which essentially covers the range 20 000 K < T < 100 000 K and n < 10^^ cm~'^. As shown in 
[4T| at high density the Pauli-Fock contributions are essential for the limitation of the number of 
energy levels and determine the transition from full ionization at very high densities to partial 
ionization. For a detailed discussion of this problem see, e.g., [T3J[5B]. 
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